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ABSTRACT 


Rainfall data from two dense rain gage networks, 
several small groups of rain gages, and National Weather 
Service (NWS) rain gages are used to describe the occurrence 
and depth-area distribution of rainfall and to develop a 
rainfall simulation program (model) for Arizona and New 
Mexico* Depth-area rainfall distribution is described b/ an 
assumed cellular thunderstorm structure, with verification 
limited to comparison of simulated and actual total rainfall 
patterns. Simulated rainfall is convenient to use in hydro- 
logic models when long term rainfall records are unavail- 
able. If point records are available, they can be used to 
verify the model* 

The program output includes accumulated seasonal rain- 
fall for any designated point (gage), point totals for indi- 
vidual events for isohyetal mapping, starting and ending 
times for all events, and Thiessen weighted watershed aver- 
ages. The output can be used directly to estimate peaks and 
volumes of runoff for very small watersheds (up to 100 hec- 
tares), and indirectly for larger watersheds with appropri- 
ate routing methods. Since the rainfall is distributed both 
in time and space, simulations of several years of record 
could be used to provide probabilities of wet and dry se- 
quences to evaluate the chances of success for range renova- 
tion programs and could aid ranchers in overall planning of 
range management programs. 

KEYWORDS: Rainfall, hydrology, model, semiarid, 
simulation, precipitation, thunder- 
storm, rain gage, range, watershed, 
computer. 
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MODEL OF TIME AND SPACE DISTRIBUTION OF RAINFALL 
IN ARIZONA AND NEW MEXICO 

By H. B. Osborn, E. D. Shirley, D. R. Davis, and R. B. Koehler 


INTRODUCTION 


Rainfall is an input used in roost hydrologic models that estiroate runoff 
volumes or peak discharge. Simulated rainfall is convenient to use in the 
development of such models because long and complete precipitation records, 
typical of areas where the model is to be applied, can be quickly generated. 
Simulated rainfall is necessary in applying such models to areas where rain 
gages are scarce or nonexistent or where historical records are too short. 

Where one or more point records are available, they can be used to verify models 
or to manipulate parameters in the models. 

In this paper, rainfall data from two dense rain gage networks, several 
small groups of rain gages (fig. 1), and National Weather Service (NWS) rain 
gages are used to develop a rainfall simulation program for Arizona, New Mexico, 
and similar semlarid regions. We first describe the rainfall occurring in these 
areas and then present a simulation program involving several rainfall models. 
Hydrologic assumptions and simplifications in the program subroutines are dis- 
cussed, and sample program output is described. 


SELECTED RELATED RAINFALL STUDIES 


Osborn and Lane (1969) studied the relative sensitivity of rainfall vari- 
ables and watershed characteristics on runoff from intense, short-duration 
thunderstorm rains. They found that, for four very small watersheds (4 ha and 
less), runoff volume was most strongly correlated to total rainfall and peak 
discharge was best correlated to maximum 15-roin rainfall. Total runoff was also 
highly correlated to maximum 15-roin rainfall, because the maximum 15-roin ra In- 
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Water Resources, University of Arizona, Tucson; and hydrologic aid, SEA-AR, 
Tucson, Ariz. 85705. 

The year in italic, when it follows the author's name, refers to Litera- 
ture Cited, p. 21 . 


1 




• Subregional Bose Stations 
0 Experimental Watersheds 

Subregional Model Borders 
Rivers 


Figure 1.— Region and subregions covered 


by regional rainfall model. 


Differing watershed characteri^lcrdid not^^dd^^ relationships, 
peak, or .ol™., runoff, hut thn u.torshedn atudLd'^erTqulu'sLIlar!' 

Fogel (i969) reported on the efferrc ^ 

-f from small watersheds in the Southwest He no/r^f variability on run 
^plicated process and becomes more c^nUea?!/ J ^ 

nse, short-duration thunderstorm ralnfrn ^ highly In- 

it methods for estimating runoff voluLs reaulr'”^W^'^ extent. Most cur 

1 depth, which can lead to total ral 

<-rors in estimating runoff. 

Duckstein et al. (1972) 

fall for summer-type storms In the so model of runoff -producing 

^hat «odlffoahlo„rf„ ru^H'L^'^ooo ^ay 

and fhnh f„ nfrhnr can. rafnfnU iaP^ruu^rh^^^opIrr^^SlX" d^l 


2 



tsrniinG actual changes in runoff with changes in watershed characteristics* 

They consider summer rainfall as an intermittent stochastic phenomenon and ob“ 
tained the probability distribution of areal rainfall by convoluting a Poisson 
distribution for the number of events with a geometric or negative binomial pro- 
bability of rainfall amount* They then used their rainfall model in several 
rainfall— runoff relationships to Illustrate the practical value of the method* 

Osborn and Laursen (1973) identified the differences in the runoff-pro- 
ducing characteristics of the airmass thunderstorms common to Arizona and west- 
ern New Mexico and the frontal convective rains more common to eastern New 
Mexico and western Texas* They found that pealt discharge from semiarid range 
land wastersheds larger than 2 km^ in southeastern Arizona was best correlated 
to maximum 30— min rainfall (as opposed to maximum 15— min rainfall for very small 
watersheds), which they referred to as the core of runoff -producing thunderstorm 
rainfall. The simulation program presented in this paper is based on the con- 
cept of a core of runoff-producing rainfall and can be used for watersheds of 
2 km^ and larger. 


RAINFALL DEPTH- AREA DURATION 


The USDA Southwest Rangeland Watershed Research Center (SRWRC) has operated 
the 150-km^ Walnut Gulch Experimental Watershed in southeastern Arizona and the 
174-km^ Alamogordo Creek Experimental Watershed in eastern New Mexico since 1954 
(fig* 1). Currently, there are about 95 and 65 recording rain gages, respec- 
tively, on the two watersheds, although not all of the gages have been operated 
during the full period of record* Also, the SRWRC operated several small exper- 
imental watersheds near Safford, Ariz., and .Albuquerque, N* Mex* , through 1976, 
and near Ft* Stanton, N. Mex., from 1966 to the present* Data collected from 
these networks are used to illustrate the range and variability of runoff-pro- 
ducing characteristics of thunderstorm rainfall. 

Regional differences in rainfall amounts and intensities in arid and semi- 
arid regions have been widely investigated; however, quantitative descriptions 
of these differences, usually as depth-duration frequencies, generally have 
ignored differences in the storm system that generated the rainfall, and have 
lumped together essentially different storm populations* For example, thunder- 
storm rainfall occurring in the arid and semiarid regions of the Southwest can 
be divided roughly into rains occurring from airmass thunderstorms and from 
frontal— convective thunderstorms. 

Differences in runoff -producing characteristics between and within the two 
thunderstorm types — airmass and frontal-convective — are best illustrated using 
data from dense rain gage networks on the Walnut Gulch and Alamogordo Creek 
watersheds. 


Airmass Thunderstorms 

The extreme spatial variability of airmass thunderstorms, which is typical 
in southeastern Arizona, is Illustrated by an Isohyetal map of total storm rain- 
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Table I. 


■Ratio of storm rainfall depth for 
far 37 storms of greater than 25 


various durations to 1-hr rainfall 
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Duration (min) 



Average 

Range 


»■« 0.76 0.88 

0.15-0.54 0.28-0.77 0.38-0.92 0.47-0.96 0.66-1.0 

sept. 10. 1967t 0.17 0.31 o .41 o.53 


1.0 


iMaxiroum 1-hr point rainfall (88 mm): the onlv j . 

greater than 75 mm on Walnut Gulch ( 1955 - 77 ), ^"“hr amount 


Depth-area curves for maximum 1-hr rainfall fr^r- »•»,„ .. 

1967, July 22, 1964, and August 17, 1957, L Walnut f September 10. 

portion of fig, 3 ) again Illustrate the limited arpal ♦- ^ ^®^®d lines in lower 
storm rainfall. areal extent of airroass thunder- 
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Frontal Convective Thunderstorms 


Frontal-convective thunderstorms, which are common in eastern New Mexico, 
can cover more area, last longer, and produce greater amounts of rainfall for 
similar durations than airraass thunderstorms. The more massive nature of some 
frontal-convective thunderstorms is illustrated ty an isohyetal map of total 
storm rainfall for the event with the maximum recorded 1-hr rainfall on Alamo- 
gordo Creek, (fig. 4). Rainfall ranged from 4 mm to 103 mm, with four gages 
recording 100 mm of rain or more. Almost all rain fell in 1 hr. Although the 
storm was more massive than an airmass thunderstorm, it was reasonably well 
centered on the watershed, and apparently extended for only short distances out- 
side the watershed. Unfortunately, the network was not large enough to encom- 
pass the entire storm. 

Ratios of rainfall depths for short durations (up to 1 hr) for 37 airraass 
and frontal-convective events of over 25 mm on Alamogordo Creek are shown in 
table 2. The average ratios of 15-rain and 30-min depths to 1-hr depths are 0.56 
and 0.81, respectively, for the 37 storms (as compared with 0.66 and 0.88 for 
Walnut Gulch). The two events with maximum 1-hr point rainfall are also shown. 
Rainfall for durations of 5 to 30 rain was about average for the June 5, 1960, 
storm and considerably below average for the August 21, 1966, storm, indicating 
that intense rainfall may last longer on some occasions than has been experi- 
enced on Walnut Gulch. Again, not enough data are available to be certain of 
this or to assign probabilities to events of the magnitude of June 5, 1960, on 
Alamogordo Creek, 


Table 2. — Ratio of storm rainfall depth for various durations to l-br 
rainfall for 37 storms of greater than 25 mm for Alamogordo Creek 





Duration 

(min) 



5 

10 

15 

20 

30 

60 

Average 

0.25 

0.41 

0.56 

0.68 

0.81 

1.0 

Range 

0.13-0.47 

0.25-0.67 

0.35-0.98 

0.45-0.99 

0.63-0.99 

— 

June 5, 19601 

0.20 

0.39 

0.54 

0.69 

0.84 

1.0 

August 21, 1966^ 

0.20 

0.34 

0.46 

0.55 

0.68 

1.0 


^Maximum 1-hr point rainfall (103 mm) on Alamogordo Creek (1955-77). 

^Maximum 1-hr point rainfall (91 mm) for second greatest event on Alamo- 
gordo Creek, 
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Figure A.—MMimum 1-hr rainfall, June 5, 1960, Alamogordo Creek. 
Each contour represents 5 mm of rainfall. 


Differences on rainfall volumes for observed maximum 1-hr rainfall are il- 

C^eerftablf maximum events on Walnut Gulch and Alamogordo 

Creek (tables 3 and 4). The more massive nature of frontal-convective storms on 
Alamogordo Creek is apparent. For example, the areal extent- In ^ 

hyets for the Walnut Gulch and Alamogordo Creek storms ^ 

12 and 148 knf . The total volume of rainfall for ^ ^ 

at least tbrpo tt,r«o c 7. „ , the Alamogordo Creek storm was 

.'iL: >» -10 only 
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Table 3 . — Maximum l~hr rainfall volumes within selected isohyets for 
storm of September 10, 1967, on Walnut Gulch 


Isohyet 
(mm ) 

Area 
(kra^ ) 

Volume 

(10® m3) 

Isohyet 

(mm) 

Area 

(km2) 

Volume 
(10^ ro^) 

80 

0.2 

0.2 

35 

38.5 

18.8 

75 

«6 

.5 

30 

47.7 

21.6 

70 

1.4 

1.1 

25 

57.0 

24.1 

65 

2.6 

1.9 

20 

^6.8 

^26.3 

60 

4.1 

2.8 

15 

^79.3 

128.5 

55 

6.7 

4.3 

10 

^ 104.0 

I 3 I .6 

50 

12.4 

7.3 

5 

1 124.0 
^ 155.0 

1 33.1 

45 

40 

22.8 

31.1 

12.2 

15.8 

0 

I33.9 


^Partial storm areas and volumes recorded only within the rain gage 
network. 


Table 4 . — Maximum 1-hr rainfall volumes within selected isohyets for storm 

of June 5, 1960, on Alamogordo Creek 


Isohyet Area Volume Isohyet Area Volume 

(mm) (km^) (106 ro3) (mm) (km^) (10 m ) 


95 

2.3 

2.2 

90 

8.0 

7.1 

85 

16.3 

14.8 

80 

25.4 

22.2 

75 

34.2 

29.1 

70 

43.8 

36.1 

65 

66.8 

51,8 

60 

97.6 

71.0 

55 

124 

86.0 

50 

148 

98.5 


45 

ll53 

llOl 

40 

ll54 

I 1 O 2 

35 

1155 

I 102 

30 

ll56 

I 102 

25 

1157 

ll03 

20 

ll59 

ll03 

15 

I 16 I 

ll03 

10 

ll62 

ll03 

5 

ll64 

ll04 

0 

ll74 

ll04 


3 Partial storm areas and volumes recorded only within the rain gage 
network. 


Depth-area curves for maximum 1-hr rainfall for the storms of June 5, 1960 
and June 16, 1966, and for maximum 6-hr storms for August 21, 1566, on Alamo- 
gordo Creek are plotted along with depth-area curves for "record" Walnut Gulch 
storms (fig. 3). The curves for Alamogordo Creek are much flatter than those 
for Walnut Gulch. Also, the longer durations plot above the shorter durations, 
indicating that, for frontal-convective events, storm areas generally increase 
with Increasing depths. 
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TIME DISTRIBUTION 


Because of the nature of alrroass thunderstorms, rainfall time distribution 
both at a point and over a small watershed (up to 150 km^) for major events has 
relatively little effect on flood peaks and volumes (Osborn and Laursen, 1973 ), 
This is primarily true because the Intense "core" of runoff -producing rainfall 
occurs as one unit and is much greater than the light rains that often follow 
and sometimes precede heavy rainfall. Based on Walnut Gulch data, intensities 
for major airmass thunderstorms are on the order of 125 mm/hr for about 30 min 
or a little longer, whereas Intensities during the mature and dissipating 
stages are generally more on the order of 2 to 3 mm/hr. Also, intense rainfall 
tends to begin and end very abruptly. 

As an example, histograms of maximum point rainfall for three of the four 
largest runoff -producing airmass thunderstorms on Walnut Gulch are shown in 
figure 5. In all four cases, including the one not shown, over 90 percent of 
the rain fell in less than 60 min, and intensities outside the core of runoff- 
producing rainfall were too low to cause runoff. In almost all cases on Walnut 
Gulch, runoff-producing point rainfall lasted less than 60 min. 


On the other hand, in eastern New Mexico, where frontal activity in combi- 
nat ion with convective heating often affects rainfall, the time distribution of 
rainfall can have a significant effect on flood peaks and volumes. Histograms 
of three major runoff-producing events are shown in figure 6. Rainfall is con- 
centrated within 60 min during two of the three events. The third event (on 
August 21, 1966) was distinctive in that Intense thunderstorm rainfall was 
superimposed on a heavy frontal-type rain with heavy runoff resulting. 


STORM MOVEMENT 


can T semiarid regions are generally ephemeral and 

can abstract large volumes of runoff in the normally dry channels. These 

relatr*^^M Complexity of thunderstorm rainfall-runoff 

effect on mH^r * ®torm movement appears to have little 

zateJLpH small watersheds. If the storm moves up the 

watershed, runoff moves more rapidly through the wetted channels, whereas if 

ar^iv Tt S'L w runoff TanL to 

at the outlet at the same time as runoff closer to the outlet Tn Mth 

E£~Hf s =r;3.-:”r“ 


REGIONAL RAINFALL MODEL 


The regional rainfall 
computer rainfall models. 


simulation model Is based 
Thp g stochastic 


primarily on two existing 
depth-area model of airmass 
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Figure 6 •“—Histograms of three major runoff— producing events on 

Alamogordo Creek. 


thunderstorm rainfall that was developed with recording rain gage records from 
the Walnut Gulch Experimental Watershed (Osborn and others, 1972 \ Osborn and 
others, 1974). The second simulates daily rainfall occurrence in Arizona and 
New Mexico (Osborn and Davis, 1977). In addition to revising and adopting these 
two models, routines were developed to simulate frontal rainfall and to expand 
t e airmass thunderstorm rainfall model to include frontal convective events. 

The regional model is a set of easily modified subroutines, so changes can be 
made in the model as more is learned about rainfall in the Southwest. 


Input Data 


Comparison of NWS rain gage records and topographic features in Arizona and 
New Mexico suggested the region could be divided into three relatively homogenous 
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subregions (fig. 1; Osborn and Davis, 2977). Subregion 1 is eastern New Mexico; 
subregion 2 is western New Mexico and northeastern Arizona, and subregion 3 is 
the remainder of Arizona. The program overestimates rainfall occurrences for 
watersheds below 300 or above 2300 m in elevation; fixed coefficients in the 
program must be adjusted for watersheds outside these elevation limits. The 
appropriate subregion must be entered along with the parameters—latitude, lon- 
gitude, and elevation — representative of the watershed. The relative coordi- 
nates of each rain gage (or dummy measuring point), along with Thlessen weights 
(optional) are also entered. 

Probabilities for airmass thunderstorm rainfall are based on seasonal prob- 
ability curves at base stations in the subregions—Roswell (base station of 
subregion 1), Las Cruces (base station of subregion 2), and Douglas (base 
station of subregion 3). The seasonal curves are adjusted automatically when 
the watershed latitude, longitude, elevation, and subregions are entered. 

Thunderstorm rainfall probabilities for SE moisture (explained in the fol- 
lowing section) are based on a 2— km^ watershed and are adjusted upward for 
larger waterseds. For example, there are, on the average, about twice as many 
thunderstorms on the 150-km^ Walnut Gulch watershed than on a 2-kro^ subwatershed 
within Walnut Gulch. 


Rainfall Occurrence 


Available recording rainfall records from USDA experimental watersheds and 
the NWS rain gage network in Arizona and New Mexico were used to develop the 
rainfall occurrence routine in the regional model (fig. 7; Osborn and Davis, 
2977). Two types of precipitation-producing systems were identified — frontal 
and airmass — with airmass divided as to the source of moisture. Moist air from 
the Gulf of Mexico is called SE moisture; moisture from the Pacific and the 
Gulf of California is called SW moisture. Frontal activity is more likely in 



Figure 7. — Simplified schematic of summer rainfall occurrence in 

Arizona and New Mexico. 
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eastern New Mexico and northern Arizona than in southern Arizona and southwest- 
ern New Mexico, and has a strong latitudinal component, whereas SE moisture is 
more common in New Mexico, and SW moisture is more common in Arizona. Both have 
a strong longitudinal component. Elevation is the third significant parameter 
in the occurrence routine. Daily rainfall occurrences are determined independ- 
ently for the three systems, and are combined, when necessary, in the depth-area 
rainfall simulation (to be explained later). 

We determined frontal rainfall occurrence by using uniform seasonal distri- 
butions based on watershed location and elevation. An element in the frontal 
subroutine increases rainfall probability on the day after it rains. The chance 
of rain is reduced to the original fraction after a dry day. The subroutine 
adjusts for frontal activity that is more common as we move both north and east 
in the region. 


We determined SW rainfall occurrence by using an annual uniform distribution 
location and elevation. With satellite pictures now available, we may be able to 
Identify significantly different seasonal probabilities in SW rainfall occurrence 
and to incorporate these differences into the program. The likelihood of SW 
rainfall occurrence is much greater in Arizona than New Mexico, and the subrou- 
tine Incorporates this difference. Also, once SW rainfall has occurred, the 
chance of rain on the following day increases (Osborn and Davis, 1977). 


The SE rainfall occurrence routine is more complex. SE moisture is sea- 
sonal, representing the so-called monsoon in the southwest. We chose base 
stations that appeared representative of each of the tlrree subregions. Daily 
probabilities are automatically modified ty a multiplier based on the watershed 
location and elevation. In general, SE moisture increases with elevation and 
decreases as we move north and west, and the subroutine allows for this. Al- 
though separate equations were developed in the model to account for persistence 
n rental activity and moisture from the southwest, no persistence equation was 
included for moisture from the southeast. 


Several investigators have used Markov chain models to predict point rain- 
fall occurrenc.. S.lth aad Schrelber (1973) assumed all avents ware of tha soma 
fitted daily rainfall occurrences at three stations In southeast- 
Arizona with a segmented first-order Markov chain model. Woolhlser (1975) 
has proposed a three-parameter mixed exponential Markov chain model of dally 
rainfall based primarily on data from the Great Plains area. Possibly this 

morp could be adapted to the southwest as a substitute for the 

more cumbersome empirical equations that are presented here. 

Moisture particularly from the southeast, dominates summer rainfall in 
pppp Arizona. Osborn and others (1972) modeled storm occurrence as a 

seasonal Bernoulli random variable based on occurrence of storms of more than 5 
mm on Walnut Gulch. Since there was no statistical difference in the persist- 
ence pattern for major events between simulated and actual data, the model did 
not Include a persistence equation for southeast moisture. 

dav thunderstorm day is simulated, more than one event may occur on that 

day. Gances for multiple events are based on Walnut Gulch data. Dates and 
ypes of rainfall occurrences, simultaneous occurrences, and multiple events 
are stored for use in the depth-area program routine. 
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Time Distribution 


Thunderstorms have an extremely complex cellular structure in time and 
space j which Is extremely difficult to model without major slmplif ications n 
Our program simulates a simplified cellular structure as a vehicle to arrive at 
storm output, which in turn can be compared directly with real storm data* 

If a thunderstorm is indicated, the program first generates a beginning 
time based on a normal distribution centered at 1700 hr* Each thunderstorm has 
between 3 and 12 cells with a mean cell distribution of 6* The cells are 
multiplied by a constant (assumed average cell rainfall duration) to determine 
the end time for the “runoff ^producing’' portion of the event (light rain may 
fall for some time after a thunderstorm, but this nonrunoff-producing rainfall 
is not included in the model)* Individual cells would not necessarily produce 
runoff, but are assumed to be in the active part of the thunderstorm event. 

If the occurrence routine calls for a second event, the program will 
repeat the routine, but only after determining that the second event begins 
after the first event ends* Beginning times are generated until this occurs* 

Depth-Area Rainfall Simulation 


The simplest output in the program occurs on a day with only frontal rain. 
Frontal rain is generated for a point and assumed uniform over the entire 
watershed# For larger watersheds, we may add a subroutine to vary rainfall 
with elevation* Frontal storms must last for at least 3 hr and can extend into 
the following day. 

Many of our assumptions for the dimensions and distribution of the inter- 
nal components of thunderstorm rainfall have been determined ty trial and error 
until the simulated and real storm outputs are comparable (Osborn and others, 
197 3) m At this time, we assume depth-area rainfall simulations for SE, SW, and 
SE/SW occurrences are the same# We assume a constant radius for alrmass thun- 
derstorm cells# In the future, we hope to be able to Identify significant dif- 
ferences in SE and SW events and to program these differences. 

The program will generate a rectangular area enclosing the rain gages (or 
dummy points) entered into the program and establish a buffer area around the 
rectangular watershed (fig* 8). The center of the first cell is randomly lo- 
cated somewhere within the buffer (which includes the watershed area). We as- 
sume the second cell is up to a cell diameter from the first cell and has a 
random directional component- The third, fourth, fifth, etc*, cells have a di 
rectional component based on the direction between the previous two cells. 

Next, cell center depths are generated independently for each cell* An example 
of the program at this point is shown in figure 9* 

Once the* center depths are determined, the rainfall is distributed and ac- 
cumulated at all designated rain gage or dummy measuring points* Average 
watershed rainfall is determined from the Thies sen weighting, and can be 
printed out in the final program# Point event totals are printed out and can 
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be used to plot Isohyetal maps (fig. 10). 

When frontal and thunderstorm events are simulated on the same day, the 
program combines the elements of both storm types. First, frontal rainfall is 
simulated and distributed to each gage. Then, the cells are simulated as de- 
scribed previously. Based on data from the Alamogordo Creek rain gage network, 
we assume that cell diameters and depths are greater for the frontal-convective 
simulation. 



Figure 8. — Rectangular watershed and buffer area created in the 
depth-area rainfall subroutine. 


EXAMPLE OF OUTPUT 


The user output can be best described by simulating rainfall in an actual 
situation. We simulated rainfall on a small (8 km^) subwatershed (W-11) on the 
Walnut Gulch Experimental Watershed near Tombstone, Arlz. 

There are 10 recording rain gages on or immediately adjacent to W-11. The 
relative location and Thiessen weight for each gage must be entered into the 
program (table 5). The program will establish a rectangular watershed and buf- 
fer (fig. 8). In practice, most gages are "dummy" locations used to identify 
the watershed boundaries and to register simulated rainfall at a point. The 
latitude, longitude, and average watershed elevation also must be entered into 
the program. The program can be instructed to simulate rainfall for an entire 
year, or any portion thereof; we have concentrated on the thunderstorm season. 

A typical season for W-11 is shown in table 6. Rainfall is concentrated in 
July and August with almost all rainfall occurring from airmass thunderstorms. 

Once it has been established that thunderstorm rainfall occurs on a partic- 
ular day based on a given probability for that day, there is a lesser chance 
of a second storm occurring, and then a third, and so on. The program will 
generate a series of rainfall cells starting within the subject watershed or 
buffer area. In our example, shown in figure 9, the program predicted two 
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X 5 rain gage 
ia=accumuloted roin 
13.0= cell center depth (mm) 


STORM I 


J 



Figure 9. — Simulated storm cell development on a Walnut Gulch subwatershed- 
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Table 6. — S imitated rainfall occurrences for a typical season for Walnut Gulch subwatershed W~21 
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^Nutaber of thunderstorms on given day. 


11. By chance, early season rainfall for both years was below average. In 
this example, the high year resulted from heavy rains in September. Based on 
5 yr of record from Walnut Gulch, such heavy September rains occur on the 
average of about once every 5 yr, and analysis of Walnut Gulch rain gage data 
indicates that September rains are not correlated with July-August rainfall. 

The program can be instructed to print out any or all designated gage accumula- 
tions as well as the areal distribution and watershed average for all events. 

The program output can be used for many purposes. For example, Isohyetal 
maps can be constructed for any event (fig. 10) or for seasonal or annual pre- 
cipitation. Thlessen weighted averages can be computed. Average or point 
watershed rainfall can be used directly as input to hydrologic models to obtain 
the probability distribution function of peak discharge and storm runoff. The 
seasonal rainfall distribution at a centrally located gage (or gages) can be 
used as input for range condition models or to estimate the probable outcome of 
range renovation programs. The latter may be particularly important, since hy- 
drologic data are usually lacking when range renovation programs are evaluated. 
Also, model output would be used, for example, to estimate the El (rainfall 
energy) factor in the Universal Soil Loss Equation (USLE) and to predict ero- 
sion and sediment yields given other appropriate models. 



Figure 11. Smoothed mean and maximum and minnimum accumulated 
seasonal rainfall for 10 yr of simulated rainfall. 


SUMMARY AND CONCLUSIONS 


Precipitation in Arizona and New Mexico occurs primarily from frontal ac- 
tivity or convective systems or a combination of both. Summer rainfall in the 
form of thunderstorms is extremely variable both in time and space, whereas 
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frontal rainfall is more evenly distributed (Osborn et al., 1979). Models of 
both the occurrence and areal extent of rainfall are necessarily simplified. 

A model has been developed that simulates rainfall occurrence and amounts 
for ungaged watersheds up to 150 km^ and elevations between 300 and 2300 m in 
Arizona and New Mexico. The model is a combination of many subroutines with a 
number of alternative Inputs and outputs. The model simulation includes accu- 
mulated seasonal rainfall for any designated rain gage, point totals for indi- 
vidual events, starting and ending times for all events, and Thiessen weighted 
averages. The model output can be used, for example, to simulate thunderstorm 
rainfall, which in turn can be used to predict the distribution and magnitude 
of peak discharges and runoff volumes for ungaged watersheds. The model output 
can also be used to estimate erosion and sediment yield from small watersheds 
and as input to more complex range renovation or management models. 

Many range management or renovation treatments proposed by ranchers rely 
heavily on rainfall as the "make or break" variable of the treatment. These 
treatments could be tested before implementation by simulating rainfall for 
varying conditions and then determining the corresponding Impact on the range 
management treatment. This preliminary testing could reduce the risk of 
treatment failure and economic and resource loss. 

The program is available by writing to; USDA Southwest Rangeland Water- 
shed Research Center, 442 East Seventh Street, Tucson, Ariz. 85705. 
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APPENDIX 


The main program, SATDOR (Space and Time Distribution of Rainfall) (fig. 
12), is broken down into a number of subroutines (fig. 13). This allows the 
operator to easily modify specific procedures simply by adding the parts of 
Interest. The function of each subroutine is explained so that the operator 
will understand its role in SATDOR. To insure consistency of the sequence, 
when the program first starts, the day before the start of the season is as- 
sumed to have no rain. The program proceeds on a consecutive daily basis. 

There are seven major subroutines; tliree are made up of minor subroutines 
(fig. 13). The first three are required to generate rainfall occurrence; the 
next two are used to generate the spatial distribution of rainfall; and the 
last two are used to determine the temporal distribution of rainfall events. 
The seven subroutines are: 

(1) WSHED - Definition of watershed boundaries and gage locations. 

(2) PDF - Probability distributions and parameters for rainfall 

occurrence defined. 

(3) DROM - (Dally Rainfall Occurrence Model) - Generates daily rain- 

fall occurrence record. 

(4) SPAR - (Space Parameters) -Spatial parameters defined and 

calculated. 
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Figure 12. — Flow chart and 
subroutine locations in 
SATDOR (Space and Time 
Distribution of Rainfall). 
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SATDOR 
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DIST GGNOF GGPOSH 

Figure 13* — Subroutine linkage In SATDOR. 


(5) SDOR - (Space Distribution of Rainfall) - Spatially distributes 

rainfall for individual events* 

(6) TPAR - (Time Parameters) - Temporal parameters defined and calculated* 

(7) TDOR - (Time Distribution of Rainfall) - Distributes rainfall 
temporally. 

WSHED 

The WSHED subroutine reads in the rain gage numbers > whether the gages are 
on the physical watershed, the x and y coordinates of each gage, and the Thiessen 
weight of each rain gage. Then, the maximum and minimum coordinates are used to 
establish a rectangluar rainfall area. Even though the physical watershed is 
usually some irregular shape, the rain gages will be used to define a rectangle. 
The rectangular area is computed and printed. A buffer of 4 km is added to each 
side of the rectangle so that simulated storms may form both outside and inside 
the watershed (fig. 8). At present, the buffer is constant, but we are investi- 
gating whether it should vary with watershed size. 
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The areal adjustment factor (AAF) Is also defined in this subroutine* The 
AAF affects the probability distribution (P) read in by PDF (see PDF section)* 
The AAF is a multiplication factor that increases rainfall probability (p) with 
Increasing area, but not in direct proportion to areal differences (fig* 14). 

AAF equals one for watersheds up to 2 ktn^, and is varied upward as rain is sijuu^ 
lated over larger watersheds. The rainfall probability cannot exceed one. 



Figure 14. — Method for adjusting SE thunderstorm rainfall 
probabilities according to watershed area. 


PDF 


The PDF subroutine reads in data used by the DROM subprograms. These data 
Include the watershed location, elevation, region, season beginning, season end 
(Julian dates), number of seasons to simulate thunderstorms, probabilities of 
SE thunderstorms for each region (adjusted hy the AAF), persistence of thunder- 
storms from one day to the next, and chances for more than one event on the 
same day. 


DROM 


The DROM subroutina contains three subroutines of its own# These are_SE, 
SW , and FRONTAL (fig. 13), which represent sources of moisture from the princi- 
ple sources, the Gulf of Mexico, the Pacific Ocean, and weather fronts, respec- 
tively. The SE and SW routines simulate thunderstorms, whereas FRONTAL simu 
lates constant low-intensity rainfall over the specified area. These systems 
can occur in any combination# Watershed location will dictate which combina- 
tions are favored# The SE and SW .thunderstorms are more massive when combined 
with frontal events. 


The SE subroutine determines whether any storms Involving this system will 
occur based on the thunderstorm probability curve for the watershed reg on, 
latitude, longitude, and elevation. The SE rainfall protebility curves a 
sites other than the base site of the region for which the ptoba es we e 
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developed are calculated based on latitude, longitude, and elevation. Also, 
the SE computations are the only parts of DROf^ that are affected by the AAF 
term (see above). 

The SVI and FRONTAL subroutines are very similar. Both have a constant 
probability for rain on a day when there is no rain on the previous day. If 
there is rain, the following day has an increased chance of rain. The chance 
of rain after two successive rainy days is decreased and so on. The equations 
used to calculate the daily probability for SW and FRONTAL have different 
constants; FRONTAL does not include watershed elevation, whereas SW does. Also, 
a constant depth is generated when FRONTAL is called, whereas SW calls for 
thunderstorm rainfall. 

SPAR 


SPAR reads in spatial parameters (cell diameter, mode distance between 
cells, number of cells) and calculates shape parameters to be used in the IMSL 
canned routines GGPOSH (a pseudorandom Poisson distributed deviate), GGAMA (a 
pseudorandom gamma distributed deviate), and GGNOF (a pseudorandom normally 
distributed deviate). Parameters and statistics are read in and printed out. 
These are used in the SDOR subroutine. 

SDOR 


The SDOR subroutine is the most complex part of the SATDOR program. There 
are four main subroutines in SDOR with several minor subroutines (fig. 13). 
Basically, SDOR generates and distributes rainfall to gages on the watershed, 
depending upon whether tlie event is a thunderstorm, frontal storm, or both. 

This is done by using a combination of subprograms particular to a storm type. 
The four main subroutines are: 

(1) DTHUND “ Generates and distributes thunderstorm rainfall. 

(2) DFRONT Generates frontal rainfall. 

(3) FDEPTH - Distributes frontal rainfall. 

(^) PRINTER Prints results of the three previous subroutines. 

DTHUND has three subroutines and one function routine of its own; DIRDIST 
and CDEPTH generate location of storm cell centers as well as depths at the 
cell centers. DIST is a function of DIRDIST and calculates distance between 
cells using GGNOF. The SPAR subroutine reads in average cell radius for dif- 
ferent storm types. The distance between cell centers cannot exceed the cell 
diameter. The number of cells per storm is generated by NOGS using GGPOSH. If 
the number is less than 3 or greater than 12, NOGS regenerates a number until 
it is within this range. At present, storm duration is a constant (for each 
storm combination) times the number of cells. GGDEPTH is a function that dis- 
tributes thunderstorm rainfall to each gage according to the distance between 
the gage and the cell centers. 

^ ppAMA frontal storm depths, and FDEPTH distributes 

no storm pattern or variation in 
errors for larger watersheds, but probr- 
jLOutine could be modified to allow for 
t, and any other topographical influence. 
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PRINTER prints results of the SDOR subroutines. 


TPAR 


TPAR reads in the temporal parameters and calculates shape parameters for 
the canned subroutine GGNOF and prints out the parameters and statistics. 
These values are used in the TDOR subroutine. 

TDOR 


Thunderstorm rainfall events are temporally distributed by TDOR using the 
functions BTTH and ETTH. The subroutine TDOR also handles the temporal distri- 
bution of frontal events. The set of statistics for the starting time, devia- 
tions, and duration will vary on storm type and if there is more than one storm 
per day. 

BTTH varies with the number of storms per day and with which storm is as- 
signed a beginning time (the second storm must occur after the first is fin- 
ished). Results are checked to see if they are valid and, if not, regenerated. 
ETTH generates the end time of thunderstorms by indirectly using 6GP0SH and the 
number of cells. 
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